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SYfTOPSIS 

TRAJ'TSITIIT □FFUCTu IN M FWCLOSIiD BODY 

DUE 70 THEDMAL DISCHARGE 

JUNE, 1979 

The study of storage of thernal energy xn 
enclosed Ajater bocb.es and Uie prediction of the perfor- 
mance under various conditions of practical interest 
ere important aspects in the design of such systems. 

The aim of the present ■'’ork to study tiie importmt 
thermal effects of a heated discharge into an endossd 
vater body. Interest, obviously, lies in dotcremining 
the rernating tei^nerature field, uirticularly t^ic groi'jth 
and decay of thermal stratification m the uatir body. 

T’lis study is also of relevance in heat rejection 
S3'‘steiis, 

In order to stud3'- the thernal effects of a 
heated d.ischargo into a ifatcr bodj^', an oxporimental arrange- 
ment uas designed and fabr.icatod, Tlio inflow rate ’^as kept 
erual to the outflo"',; so as to maintain a stead.;>' flo'i; situa- 
tion. For given inlet conditions, tlie trensient behaviour 
of the ^wfa/cor, in tens of tli? outflou tciiperature and the 
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temperature profiles measured in the tank, was studied at 
various locations, by means of thermocouples. Transient 
temperature profiles were obtained during the discharge of 
hot water. The decay in the temperature profiles was also 
studied in the absence of flow. It was observed that the 
transient response of the enclosed water boc3y is strongly 
dependent on the configuration of flov/ circuit. The growth 
of thermal stratification in the tank with time was studied 
and the observed temperature profiles indicated strong 
mixing in the region above the inflow. 

An analysis has also been carried out to determine 
the nature of the flow and temperature fields in the tank. 
Various mathematical models were formulated in order to 
evaluate the transient response of the enclosed vjater body, 
under various conditions relevant to practical sj’-stems and 
to the experimental arraiigement, 

A simplified zero dimensional model, with uniform 
temperature in the tank, was initially employed to predict 
the performance of the water body under various heat trans- 
fer conditions at the surface of the tank. The energy 
equation yielded the effect of thermal losses on the tran- 
sient behaviour. The results obtained were found to be 
in reasonable agreement v/ith the observed experimental 
results, under conditions which gave rise to considerable 
mixing in the water body. 



A one dimensional model was then formulatedj 
assijming complete mixing in the horizontal plane. The 
governing equation which included an averaged vertical 
flow was solved by finite difference techniques. The 
results have been compared with the temperature profiles 
obtained experimentally. The model was later modified to 
incorporate the vertical mixing due to the buoyancy effects 
and the flow. The modified model was found to yield re- 
sults which were m good agreement with the experimental 
observations. 

The analysis was also extended to a two-dimensional 
protlem in order to determine tlie nature of flow field that 
arises. The cases of inviscid flow and of creeping flov/ 
were initially/- considered for various inflow and outflow 
locations. Streamlines were obtained for these various 
configurations. The flow field thus obtained indicates 
qualitatively the effectiveness of various flov/ configura- 
tions in energy dissipation as well as in thermal storage. 

For VISCOUS flow, normalized vorticity transport 
and stream-function equations were obtained. These equa- 
tions were solved by means of finite difference techniques 
for the relevant boundary conditions. The corresponding 
energy equation, which is coupled to the flow through 
buoyancy, was also obtained and solved simultaneously, TIhe 
alternating, direction implicit (ADI) method of Peaceman 
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and Badiford was en:5)loyed. The problem could not be solved 
for the actual Reynolds-numbers encountered in the experi- 
ments due to convergence problems. It was observed that the 
computer time required to solve such a coupled system of 
equations is veiy high and the approach to the steady state 
is extremely slow. The streamlines and isotherms were 
obtained. 

Several important conclusions can be drawn from the 
present work. The simplified zero dimensional model, which 
assiimes uniform temperature in the water body, can easily be 
applied for estimating the losses and in designing the system 
for a well mixed water body. Although ttie flow in a water 
body is complex and depends on the geometric locations of 
inflow and outflow, the one dimensional model can predict, 
with accuracy in many cases, the teiT5)erature profiles in a 
thermally stratified medium. The two dimensional model allows 
a determination of the flow field and may often be used in 
conjunction with sin5)lified one dimensional models for provi- 
ding information necessary for the design of such systems. 

The study would help in determining the locations of the 
inlet and discharge points for maximum effectiveness of the 
system and would assist in determining the level of thermal 
stratification that anses in an actual system. 



CHAPTER 1 


INTBDnJCTION 


1.1 WATER BODIES 

Recent years have seen a growing concern v/ith 
energy, its generation, efficient utilization and the 
ultimate disposal of waste heat. Heat rejection is a matter 
of considerable concern since it affects the environment, to 
which the energy is rejected, and often also affects the 
efficiency of the basic energy generation system. An immense 
amount of heat is rejected to water bodies by nuclear and 
thermal power plants and by chemical and manufacturing indus- 
tries. Thermal discharge into a water body might give rise 
to ecological problems. Enclosed water bodies are also used 
in the storage of thermal energy, as sensible heat, in solar 
energy systems. As a consequence, a study of the thermal 
effects arising due to thermal discharge to a water body is 
of considerable importance in heat rejection and energy 
storage systems. 

An analysis of the system involves a study of the 
flow and ten^jerature fields generated in an enclosed water 
body due to thermal discharge. This is a very coH5)lex prob- 
lem, due to the nonlinear nature of the partial differential 
equations that govern the physical mechanisms. The flow and 
the thermal stratification that arise depend on the loca- 
tions of the inflow and outflow and on various physical 
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variables, such as flow rates, ten^erature levels, dimensions 
of the water body, etc. The mathematical formulation of the 
flow and the heat transfer problem is again con^jlicated 
because of generally coupled nature of the velocity and tem- 
perature fields. They are coupled through the buoyancy 
mechanisms that arise due to teii^erature differences. The 
present investigation is directed at determining the inqpor- 
tant thermal effects of a heated discharge into an enclosed 
Water body. Interest lies in studying the ter^erature field, 
particularly the growth and decay of thermal stratification 
in the water body, 

1,2 LITEBATPBE SURVEY 

The problem of thermal discharge to a water body 
has been of interest to several investigations, largely 
because of the concern with the resulting thermal effects, 
in the water body, and, recently because of the interest in 
thermal energy storage in solar energy systems. Various 
mathematical models have been developed for the energy 
exchange mechanisms at the surface of natural water bodies, 
see, for instance, the work of Eafhael [ 1962 3, Hindley and 
Miner! 1972 3 * Henson et, al. [ 1961 3 , etc,, end for ' 
the convective, usually turbulent, processes within the 
water body, A review of the work done on heat rejection to 
natural water bodies, such as lakes, has been earned out 
by Jirka, et. al 1197^1 . 



3 


14oore and Jaluna [1971, 19721 considered the 
thermal effects of power plants on lakes, in terms of its 
temperature cycle. The stratification effects in lakes, 
during Ghe yearly cycle, were studied by Dake and Harleman 
[1969]. Spalding [ 1977 ] discusses the basic mechanisms that 
govern the behaviour of pollutant disposal in river, bays, 
lakes and estuaries. Various physical and mathematical 
models are discussed. 

Several papers on one dimensional mixing in turbu- 
lent stratified fluid flows are given in the book edited by 
Spalding and Afgan [1977 ] • The formation of an isothermal 
layer in a lake m Japan, the development and erosion of the 
thermocline, warm water spreading in cooling ponds of steam 
power plants and the development of the turbulent mixed 
region in a stratified medium have been discussed. The 
effect of turbulence cn flow in stratified fluid is also 
considered, 

Cabelli [1976 ] carried out a numerical study of 
the flow and tenqoerature fields in a thermal energy storage 
tank for various geometric configurations. The effect of 
a variation m the entrance Reynolds niaraber and the contri- 
bution of buoyancy mechanisms in promoting stratification 
have been examined. A comparison was also made with 
simple one dimensional models. The problem was solved for 
low Reynolds numbers, upto 200 (based on the perimeter of 
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the tank) , although the mixed convection parameter Gr/Re^ 
was in the range expected in practice. 

A Verification of the water temperature predictions 
for deep stratified reservoirs has boon carried out by Burt 
[ 197^^ • Numerical modelling of the thermal stratification 
in a reservoir with large discharge to volume ratio has been 
performed by Park et. al. [ 1973 3 . The predicted ten^erature 
profiles were coii 5 )ared with the field data for two different 
years. It was found that diffusion was the predominant heat 
transfer mechanism in the hypoliminion or middle layers y 
while the surface effects dominated the epliminion or upper 
layers. Oberkampt and Crw [ 197^3 have studied the velocity 
and temperature fields due to flow in a reservoir. The 
effects of inflow and outflow, wind shear and the energy loss 
on the reservoir are discussed. 

Mollondoil' and Gobhart [ 1973 A] considered the 
effect of thermal buoyancy on ^et stability. Vertical 
buoyant jots were considered and an elaborate experimental 
system was designed to study the stability of the flow. The 
temperature field was studied by means of a Sc±ilieren 
system, A perturbation analysis was also carried out by 
MoUendotf and Gebha±*t[1973 B] on the base flow to study the 
effect of buoyancy on the resulting flew. Fay [ 19731 con- 
sidered buoyant plumes in detail. laluria and Gebhart 
[ 1973 3 assessed the stability of buoyancy induced flows in 



a stratified medium. Extensive experimentation was earned 
out to obtain the desired form of stratification and to study 
its decay. Turner [1973 1 gives a review of the work done 
on similar buoyancy induced flows. Thermal stratification 
IS considered in detail. These papers, along with the seve- 
ral others on stratification and jet flows [see, Gebhart, 

19791 , indicate the interest that exists in the problem of 
thermal discharge to water bodies. Thermal discharge is 
largely by means of turbulent jet flows and the water bodies 
are generally thermally stratified. 

1 ,3 PKSSMT WOBIC 

Very little expenmental work has been done on 
the thermal transient effects of an enclosed water body duo 
to a thermal discharge. The aim of the present investigation 
was to carryout such an experimental study for large dis- 
charge to volume ratios and to verify the validity of the 
various mathematical models, in terms of the experimental 
results obtained. Two dimensional models, for the transport 
of momentum and energy in mixed convection flow, are consi- 
dered, in addition to simpler isothermal and one -dimensional 
models. 

An experimental arrangement was designed and 
fabricated to study the thermal performance of an enclosed 
water body. The inflow and outflow rates were kept the 
same, to maintain a fixed voliime of water. The effect of 
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a variation in the inflow and outflw locations on the 
thermal stratification and on the transient behaviour of the 
temperature profiles in the tank have been studied m detail. 
The results are presented in teims of the measured ten^jera- 
ture field. The buoyant flow has also been visually observed 
by means of a die injected into the tank through the inflow 
of water. 

Mathematical models have been developed and related 
to the experimental results, in order to predict the perfor- 
mance of the water body in actual systems. An isothermal 
model was employed for a fully mixed water body and its energy 
exchange with the environment has boon considered. The results 
obtained were found to be in agreement id-th the observed 
experimental results for the conditions which give nsc to 
complete mixing in the tank, particularly the case when the 
inflow is at the bottom of the tank and outflow at the top. A 
one dimensional model was also developed and predictions for 
the temperature profiles wore obtained. The model was later 
modified to incorporate the vortical mixing due to buoyancy 
effects and the flow. Very good agreement with the experi- 
mental results was obtained. 

The flow fields for the two cases of inviscid and 
creeping flow were also obtained for various configurations 
in steady flow. For viscous flow the normalized vorticity 
transport and stream funcbion equations were solved numeri- 
cally, The coupling with the energy equation was also 
considered and results vrere obtained for various values 
of the Beynoids number. 



CHAPTER 2 


EXPEKEI4ENTAL STUDY 

2,1 EXPERIMENTAL ARRAI^TGEMENT 

In order to study the thermal effects of a heated 
water discharge into an enclosed water body, an experimental 
arrangement which incorporates the basic features involved 
in the study was devised. Of particular interest in the 
present investigation was the deternination of the basic 
characteristics of tiie thermal stratification, that arises 
for various inflow and outflow conditions. The experiment 
employed an enclosed water body contained in a tank, in 
which hot water, heated by means of a constant water bath, 
was discharged. The tenperature field was measured by 
means of thermocouples, emplo 3 ia.ng a ^8 channel recorder, 
and the flow rate was obtained from suitable flow meters. 

The photographs of the main equipment used in the 
experimental arrangement are shown in Fig. 2.1. The water 
body was simulated experimentally by means of water contained 
in a perspex tank, of inside dimensions 0,6 x 0,6 x 0.^4 m, 
as shown m Fig, 2,'!. (a). On one side of this tank, three 
circular holes were drilled and brass tubes of internal 
diameter 0,6 cm, were inserted in these holes and sealed 
against leakage. On the opposite side of the tank, an 
outflow opening was made through which the outflow was 
kept equal to the inflow. A thermocole sheet, of size 
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0.6 X 0,6 X ,023 m., was used as the lid for the tank, as 
shown in Fig. 2,2 (b) . Holes were made in the lid for the 
insertion of the thermocouple probes for recording the 
temperature profiles at various locations. The desired 
inflow of hot water, at various inlet ten^jerature levels, 
was obtained by the flow of water through a coiled tube 
which was located in the constant temperature bath (ultra 
thermo-state, ten^j. range - 30 °C to 200 °C, capacity 15 
litres, control accuracy j: 0.03 ^C). The temperature of the 
hot water was varied by operating a thermoregulabor valve 
which controls the bath tenperature. In order to vaiy and 
measure the hot water flow rate, a rotameter, with a needle 
valve, was fitted at one end of the tube which carries the 
water through the bath. The flow rate was also checked by 
collecting the outlet water in a measuring flask several 
times during the experiment. The water level was kept 
unaltered, 

A Hone 3 n/ell forty-eight point temperature recorder 
was used to obtain the ternperature at any point in the water 
body. Standard copper-constantan thermocouples, made from 
wires of diameter 0,125 mm., were calibrated and used in the 
measurement of the temperature field. The recorder had an 
accuracy of the order of 0,15 ^G, A thermocouple probe as 
shown in Fig, 2,2 (c) was prepared. A scale was fixed with 
the probe and holes 1 cm, apart were drilled along its leng'ih. 
These holes were used for locating the thermocouples. The 



probe cmjld be positioned at any desired location in tlie 
tank and the flow could be varied over a \ 7 ide range of flow 
rates. The transient temperature nrofUes were measured by 
means of tins arrangement. 

The experimental arrangement employed for the flov/ 
is shown schema-cicaD ly in Fig, P.,3, Water from an overhead 
tank enters the constant teraperature bath coil at one end 
and the hot water leaving the coil was discharged into the 
enclosed vie^ter body at any one of the desired locations, a, 
b or c, as shown in the diagrair.. The outflo^r from the tank 
was taken from location d or a.ay one of the remaining 
locations, designated as a, b or c, VJhen the inflow and 
outflow locations were selected on the same side of the 
tank, the discharge from the tan]: was allowed to flo \7 through 
a flow-raoter v/ith a needle valve. Tins krrangeraent facili- 
tated the regulation of flovj, 

Thormocouplos a/ ere fitted at each inlet and ouiiLet 
location. Temperatures at other desired locations were 
obtained from the thermocouple probes descj’ibed above. The 
lnflo^7 and outflow water temperatures and the temperature at 
various other locations wero recorded at different time 
intervals. Experiments vrere conducted for various g eome- 
tric configurations, different flow rates and various 
inlet temporaturos. Extensive experimentation was also 
earned out to check the repeatability and accuracy of the 
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experiments, A very high level of accuracy, within a few 
percent, in the temperature data was obtained. The resiiLts 
obtained indicate the transient thermal behaviour of the 
water body to the energy irput, through a heated discharge. 
Some of the important and characteristic results are presen- 
ted in the next section. 

2,2 EXPEBI MENTAL RESULTS 

In the present investigation, transient teii: 5 )era— 
ture distribution in the enclosed water body has been studied 
for various flow rates, vjith particular interest in the basic 
characteristics of the thermal stratification that arises for 
different inflov; outflov^ locations. The heated, isottiemal 
or stratified, water body was also allowed to cool in the 
absence of flow and the temperature decay was studied. The 
following dimensionless variables are employed m present- 
ing the results ; 



where t is the water temperature at any point in the 
tank, t is 1±ie ambient ten 5 )erature and t^ is the inlet 
hot water ten^erature. Also, z is the vertical distance 
of the location where measurements are tal^en, from the top 
of the tank, d is the depth of water in the tank, ® is 
the thermal diffusivity of water and t is time in seconds. 
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The exit temperature is denoted by t^ (where the over 
bars on dimensionless time has been dropped boyond this,) 

The transient response of the water outlet terape~ 
rature is shovm in Fig, 2,lf with the hot water inlet at 
the bottom and the outlet located at the top of the tank on 
the opposite side. It was found that complete mixing occurs 
in the tank due to the flow and the buoyancy effects, since 
the water body was found to be isothermal for idle flow rates 
considered. The value of the mixed convection parameter is 
of the order of 0,01, based on the inlet tube diameter, 
indicating dominant forced convection mechanisms. 

The transient teii5)e rature decay of an initially 
isothemnal water body is shovm in Fig, 2.5* The average j 
loss coefficient has been obtained graphically as shown in 
Fig, 2,6, The water body remains isothermal while cooling. 
The average film coefficient h was estimated to be 16,1 
W/n?°K assuming that the top and bottom of the tank are 
perfectly insulated 

Figures 2,7, 2,8 and 2,9 indicate the transient 
theimal behaviour of the enclosed water body for the inflow 
and outflow conditions shown. The transient temperature 
response of the exit water temperature is indicated in 
Fig, 2*7, The temperature is seen to rise gradually to the 
inlet temperature level. Figure 2,8 shows the variation 
of the teirperature profiles inside the tank during thermal 
discharge. It is evident from these curves that con^lete 
mixing occurs in the tank from the intake location to the 
top of the tank. It is due to the rise of the buoyant 
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fluid and the mixing due to the flou that is generated. 
Figure 2,9 shows the decay of the teii 5 )erature profiles in 
the absence of thermal discharge, v;hen the enclosed water 
body was allowed to cool in the natural environment of the 
room. Since the cooling is largely by natural convection 
and as water has a high thermal capacity, the cooling rate 
was very slow as expected. VJhen the inlet and outlet were 
both located at the top of the tank, on opposite sides, the 
thermal behaviour of the water body was studied and is shown 
in Figs, 2.10 and 2.11, Figure 2,10 clearly indicates the 
rapid rate at which the outflo\7 temperature approaches the 
temperature of the inlet water. Figure 2.11 shows the 
thermal stratification that arises during the thermal dis- 
charge* The enclosed water body is obviously strongly 
stratified. The bottom remains largely unaltered and a 
strong teirperature gradient is generated, which becomes 
steeper with time. Figure 2.12 coiipares the transient 
response of the outlet water temperature obtained for 
various flow configurations. As expected, the inflow at 
the top gives the fastest response. 

Figures 2,13 to 2.22 indicate the thermal beha- 
viour of the enclosed water body due to heated discharge, 
when the inflow and outflow are both located on the same 
side of the water body The "Reynolds number based on inflow 
tube diameter is equal to 2^00, ’ Figure 2*13 shows the 
transient behaviour of the outlet water temperature. 
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Again the temperature is seen to rise at a decreasing rate 
to the inlet temperature value. It takes about 6 hours to 
reach v;ithin 67 percent of the inlet temperature. Figure 
2.11+ snows the teraperature profiles during thermal dis- 
charge, It is seen from these curves that the water body 
temperature at the bottom is essentially unaltered by the 
heated effluent for a fairly long time, about seven hours. 
The bottom layers gain energy by thermal diffusion and, as 
such, the rate of energy transfer is low compared to the 
top layers. The upper layers are found to be largely iso- 
thermal due to the mixing that arises. Figure 2.15 shovrs 
the decay of temperature in a stratified water body in tiie 
absence of flow. The profi3.es are very interesting. Again 
a slow decay is observed, with the \7ater body indicating 
an isothermal upper layer and a deca^^ing stratification 
level with time. The surface temperature decreases sharply 
and the bottom temperature increases due to thermal diffu- 
sion. The variation in the outflow, water temperature and 
in the ten^erature profiles with time are shown in Figs. 

2.16 and 2,17 respectively, for inflow located at the center 
(location b) of one side and the outflow at the bottom of 
the same side (location c). The outflow temperature is 
not affected by the thermal discharge for 160 minutes. 

The stratified temperature profiles are set up at a very 
fast rate, again with an isothermal upper layer, The 
bottom ten^erature nses more gradually with time. 
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Figures 2.18 to 2,22 slio\^ the effect of various 
other geometric configurations, of the inflow and outflow, 
commonly encountered in heat rejection systems. Figure 
2,18 clearly shows that, for a long time the outflow ten^je- 
rature is not affected by the heated effluent but as time 
progresses it is seen to increase. Tills observation is of 
relevance to the design of a heat rejection system, eii 5 )loy- 
ing a water body, since the outflov; is generally en 5 )loyed 
in the heat rejection process and a higher temperature of 
outflow will reduce the efficiency of the system. 

The temperature profiles at eight time levels 
are shoirm in Fig. 2,19. It is evident that the exit tempe- 
rature is not affected for five hours but is significantly 
affected as time progresses. The profiles are also strongly 
time dependent. 

For inflow located at the top of one side and the 
outflow at the center on the same side the response of the 
enclosed water body is seen in Figs, 2,20, 2,21 and 2,22, 

The variation in the exit water temperature is indicated in 
Fig. 2.20, The rise in the exit water temperature is slow 
up to 150 minutes beyond which the trend of the ciuve 
coHpletely changes. After eight hours the exit tempera- 
ture was found to be only 65 % of the highest temperature 
(i.e. inlet water ten 5 )erature ) . It is evident from Fig. 

2.21 that the bottom teii 5 )erature of the water body is not 
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siljmficr.ibl;’' affccbjc!. by ’T.a''Cjd affluenbs aven afbar 
eiglib boi-Tj, ibis inclicates the slow rate of energy 
transfer to the bottom layers. Figure 2.2 exhibits the 
cooling of the stratified enclosed water body in the absence 
of flow, when it was allowed to lose energy in the natural 
environment of the room. The bottom ter/^jerature of the 
water body rises due to thermal diffusion but at a very slow 
rate as expected. The surface temperature falls gradually. 

Photographs for the heated discharge have also 
been taken by injecting dye tracer at the inflow. Figure 
2,23 shows these photographs, as obtained from the side of 
the tank, for various inflow and outflow locations. Figure 
2,23 (a) depicts the trajectory of the discharged water 
jet when the inflow and the enclosed water body are at the 
same temperature, i.e, without buoyancy. Figures 2,23 (b), 
(c), and (d) have been taken for different inflow and out- 
flow configurations. The vertical rise of the heated fluid 
due to buoyancy and the resulting mixing is clearly seen, 

2.3 CONGLUDIMG 

For the flow under study, turbulent mixing is 
expected to have a significant effect on the energy 
exchange in the water body. In the present case, the 
buoyancy force is vertical and the thermal discharge is 
horizontal. As such, a horizontal homogeneity in the 
water body was seen to arise, 14ixing due to flow also 
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causes the region above the inlet to he isothemal. The 
flow In this region is corqplex and three dimensional. 
However, the thermal effects are found to be mainly one 
dimensional for the flow rates encountered in the experi- 
ment, The experimental study indicates tnat, for the 
conditions encountered in the experiments, it may be possi- 
ble to estimate the thermal performance of the enclosed 
water body by suitable one dimensional mathematical models. 



CHAPTER 3 


ANALYTICAL MODELS 


3.1 INTRODUCTION 


The thermal behaviour of a water bocfc^ depends on 
the shape and dimensions of the water body, configuration 
of the flow circuit, temperature and flow rate of discharge 
and energy loss to the environment. The buoyancy driven 
fluid flow and energy exchange mechanisms, which arise in 
the water body, are very complex in nature and the result- 
ing processes have, therefore, not been studied in desired 
detail so far. The relevant equations which govern the 
flow are difficult to solve without considerable simplify- 
ing assuii 5 )tions. The basic problem is a three dimensional 
one, in which the flow is coupled to the energy equation, 
through buoyancy. The partial non-linear nature of these 
equations makes an analytical solution possible only for 
very simplified circumstances, which are often physically 
unrealistic. Even the numerical solution is fairly invol- 
ved for the general problem. Therefore, some simplified 
models are developed here. Their applicability to actual 
systems is considered in terms of the experimental results 
given earlier. 
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3.2 ISOTHEBl'lAL MODELS 

In certain circiirastances, complete mixing In the 
horizontal and vertical planes of the water “body is obtained. 
In the experimental study undertaken for the present work 
also, it was observed that for inflow at the bottom of the 
tank, V7ith outflow at the top, an unstratified water body 
was obtained. Similarly for small tanks with large mixing 
effects due to buoyancy and the flow, an isothermal assimp- 
tion may be made. It is then possible to predict the 
thermal performance and to study the effect of losses on the 
ten^jerature of the v/ater body, without considering the 
detailed mechanism of the flov/. 

The equations v/hich govern the performance of an 
isothermal water body arise simply from Ihe conservation of 
energy. The following assumptions are made: 

1 , There is complete mixing m the water body and the 
volume of water remains constant. 

2, The heat capacity of the fluid is much larger con^ared 
to that of the tank material. 

3* The physical properties of the fluid are constant. 

4. The ambient temperature and the film coefficient are 
taken as constant. 
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3.2.1 Constant Inflow lomDerature 

Consider a water body as sho\m m Fig. 3*1 (a). 

Let the temperature of \:ater and that at the inside surface 

of the contemer v/all be t at time t • Water enters the 

tank at a constant temperature t^ and leaves it at the 

temperature of the water body. An onerg^^ balance of the 
water body gives : 

P C V ^ = m C„ (t. - t ) - Q. (3.1) 

p uT pi 1 

where m is the mass flow rate, Cp is the specific heat 
and P IS the density of water. Also Q-^ is the total 
heat loss from the surface. A convenient and efficient way 
for expressing Q.j_ from the outer surface of the water 
body IS given by the steady heat loss expression: 

^ ^ T'tW ^ -^s " ^a^ 

where A is surface area for heat loss. Here it has been 
s 

assumed that energy loss is proportional to the temperature 
difference (t - t^) and the Biot number EL = h dx/K 
where dx is the thickness of the wall material and Kgff is 
its effective thermal conductivity. This expression takes 
into account the heat loss by conduction from the walls of 
the water body. 

L.t T = 

The mass flow m = P i where is volume flow 


rate. 
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Inserting the above values in Equation (3.1), 


(Dt 


Q 

“1 

V 



h 

, s 
?>VQ, 


1 

1 + Bi 


) } T 


(3.2) 


Equation (3*2) is an ordinary differential equa- 
tion ■\7ith constant coofiicients axid can, therefore, be solved 
conveniently, for the initial condition T = 0atT = 0. The 
solution of Equation (3.2) can be written in the fonn 


T = CS (1 - exp~® ) 


(3.3) 


^ K W ^ 

where CS = 1 / (1 + - " •" ^ 

^1 p 


0 = ( 


h A Li, 


P V C. 
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) T , such that 


0.5< CS£ 1, and Loss Coefficient 
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1 + Bi 
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h dx 


9 S'i ;h that, 
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eff 


0 £ h 


~ jt V\ =. Vj • \»(« 


3.2,2 Constant Temnorature Difference 


In this case, the water is supplied at a constant 
temperature difference At above the outlet water ten^je** 
rature, as shown in Fig. 3,1 (b). This implies a constant 
energy input into the system, by means of, say, solar 
energy or an electrical heating system. The rest of the 
specifications are the same as before end the governing 
equation is obtained as 



i = PQ^C 
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(3,4) 


where n = h • L^, 

t - t 

I = -FT- 

Inscrbing the dimensionless variable T for t, 
Eg,. (3*4) reduces to 


dT 


!k 

V 


h A. 


PV c 


(3.5) 


For the initial condition T = 

IS 

T = CS (1 - Exp"® ) 


0 at T = 0, the solution 


(3.6) 


whe re 
1 < CS < « 


^ ^i 

CS = - - and 6 

h 


h A. 


’-pY'c — "fj 


3,2,3 Natural Cooling 

Consider a water body which is at a temperature 
% at T =0. There is no flo\^ in the tank and the water 
body IS allowed to cool in the natural environment, as 
shown in Fig, 3,1 (c). This circumstance is of relevance 
in the determination of temperature decay for energy storage 
systems, such as those en 5 )loyed in chemical and allied 
industries. The governing equation is 

0 V 0 al = - E 

t - t. 

Nondimensionalizing by T = 


equation reduces to 
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OU. _ s m 

dT "■ p V C 

P 

For the initial condition T 
solution IS 

= Exp"® 

= ^ T 

PVGp 

3*2,4 Besialts of the Isothermal Model s 

Figure 3*2 indicates the variation of ten^je nature 
T/CS v;ith time 6 • From this graph, ilie time required to 

attain a certain teirperaturo level for a given thoimal dis- 
charge can be determined. The time required to attain 99^ 
of the maximum temperature is obtained as 6 = 4,6, 

The system constant CS versus the loss coeffi- 
cient for a given inflow is shovm in Fig. 3,3 for various 
values of the heat transfer coefficient h (Q^ = 0,1 x 

o 2 

mvsec and A = 6,0 m )• It is seen that for h = 10 W/ 
s 

m^ °k » system constant reduces to 0.87, vdiereaa for 
h = 100 W/m^ *^15: » it attains the minimum value 0,5 at loss 
coefficient equal to 0,7, System constant versus the loss 
coefficient for a constant temperature difference & % is 
indicated in Fig. 3.4, for flow rate equal to 0,1 x 10*"^ 
m^sec. In this case, system constant is directly propor- 
tional to the flow rate and inversely proportional to the 
loss coefficient. 


T 

where 

9 


( 3 . 8 ) 

= 1 at T = 0, the 

( 3 . 9 ) 
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111 Fig. 3. 5, tliG transient deca*' of an isothemal 
water body is considered. It is an exponential decay. For 
any given systesj, tlio docreaso in ten'^ieraturo during a ccr~ 
tain time interval can bo estimated from bins graph. 

In Fig, 3*6, the experimentally measured ten^iera- 
turc rise in an isothermal enclosed water body is compared 
with the theoretically determined values for a constant 
inflow ten 5 )erature , A fairly good comparison betv/ecn tlie 
two IS obtained, for h = 16,1 W/m'^ °C, leading support 

to this model for systems where complete mixing arises, 

3.3 ONE Dll^SIOm IDDELS 

The experimental study indicated that, despite 
the flow in the water body being three dimensional, the 
measured ten^erature profiles that arise are essentially 
one dimensional, for the various geometric configurations 
and the flow rates considered, Cabolli f 1976 J carried out 
a numerical study of a storage tank and he also concluded 
that the thermal^ performance of the tank can be easily 
approximated by simple one dimensional models. In this 
section, suitable one dimensional models have been consi- 
dered and analyzed. 

Consider the water body as shown in Fig, 3*7 (a), 
in -viiich the inflow is located at a certain depth and 
outflow is at the top. The energy equation is obtained for 
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a vertical one-dmiensional nodel as 


3 1 , at 
sT + " 3^ = “ 


A P C. 


(t - t,) 


( 3 . 10 ) 


where w is the average vertical velocity, a the 
molecular thermal diffusivity* P the perimeter of tlie 
water body and A its cross-sectional aroa. Here loss at 
the sides is considered iii terms of h, which may as before 
be taken equal to h = h (1/1 + Bi) . In writing the expre- 
ssion for the heat loss from the sides of the water body, it 
has been assumed that the film coefficient h also takes 
into account the radiative heat loss. In general, however, 
the convective and radiative terns may be employed in full. 


The equation (3»10) can be non-dimensional ized by 


defining 


T = 


% “ ^a ’ 


Z = f , 


where d is the depth of the tanl:. 

The nondimensional equation therefore becomes 


3 T ■ ^ AT, 

TT ^ 9 z 


- HT 


(3.11) 


where the overbar" on dimensionless time has been doropped 


for convenience. Also, 


H = 


h P d 


h P 

f\ oi & C|p 
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If we assume an initial condition of uniform tempe- 
rature and a step change in the inflow water temperaturs , 
the boundary conditions are 

T ( 2, 0) = 0 

T (2^, T ) = 1 

= 0 at 2 = 0 and 2=1 

9 2 . ^,^7 

In section 11, vertically below the inlet, heat 
transfer is largely by diffusion, convective mixing being 
essentially negligible. For a one dimensional model, of 
course, no velocity exists below the inlet point, the flow 
being above the inlet to the exit at the top. The abcve 
equation, therefore, reduces- to the following,, in the region 
below the inlet, 

« HT C3.12) 

3 T ^ 2r 

The initial and the boandary conditions are 

T (2 , 0) = 0 

T (2j,, i; ) = 1 

=0 at Z = 1 2 = 0 

3 Z 

Here, it has been assumed that the to 2 _an^ 
bottom of the body of water Is perfeottytos^^d. Hie 
gradient may be obtained from energy transfer, If heat 
loss IS considered. 
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3*3*1 ! Melliod of Solution and Results 

Equations (3.11) and (3.12), with the given 
houndarj," conditions, v/ere solved numerically by finite 
difference techniques by employing an implicit method. The 
tank was divided into sections, each of uniform teirpera- 
ture. This produced a tndiagona]. system of equations which 
was solved by a tridiagonal algorithm, described by Carnahan, 
et, al. { 1969 J • 

For inflow at the mid point and the exit at the top, 

the numerically obtained temperature profiles, for the flow 

— Cj J4- 

rates of 0,67 x 10 mvsec and 0,1675 x 10 mvsec, are 
shown in Figs, 3.8 and 3*9. The predicted temperature dis- 
tribution IS clearly unstable, since a temperature increase 
with depth IS observed. This causes denser fluid to be above 
the lifter fluid and is, therefore, expected to be thermally 
ims table. In practical situations, this type of a tempera- 
ture profile will generally cause mixing. It may exist only 
when thv- losses from the surface are at a very high rate,* 
as indicated by Drake and Harlemann [ I 963 I . However, they 
pointed oub that a slight disturbance will cause vortical 
mixing. This vortical mixing vrill take place up to a finite 
depth Dj^ and can be determined by making an energy balance 
(Refer, Fig. 3.7 (b)) as 

* Author also observed this experimentally with 
thb water body uncovered at the top . 
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J { (T (Z, T ) - Tj^) }dZ = 0 (3.13) 

o 

at z = , I (D^, r ) = T_^ (3.1W 

where is the temperature of isothermally mixed layer. 

From Equations (3.13) and (3.1^), the depth of 
the unstratified mixed layer near the surface of the water 
hody, caused due to the mixing arising from the unstable 
thermal stratification, can he determined. In a buoyant 
flow, as in plumes, a Ingher temperature is expected at 
greater depth. However, in the experiments carried out here, 
buoyancy induced mixing and that due to the flow was seen to 
give rise to an isothermal upper layer, 

3,3*2 Isothermal Unner Pegion 

The experimental results showed tliat, for the 
range of flow rates studied, if the outlet is situated at 
the top and inflow at any other height, then complete mixing 
IS observed from the inflow to the top of idle tank. Thermal 
energy is transported by the convective process, ■vdiich also 
includes the buoyancy mechanisms. Below the inlet, the 
convective flow is absent and conductive mechanisms domi- 
nate. 

Consider the water body as shown in Fig . 3*7 Cc) • 
From the inflow to the outlet at the top, the water body 
is assumed to be isothermal at any instant. The overall 



energy balance for the isothermal region, for unit horizontal 
area, gives 


3 1 ^ V 

ax ~ 
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a__ ^ 

Z^ 3 2 


(t-V (3.15) 


where is the depth of the isothermal region and is 
evaluated at the interface of the t\7o regions. 


If we nondimensionalize as before, the equation 

becomes 
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( 3 . 16 ) 


where 



W = 


SLd 

a 


h P 

= AK^ 


Equation (3»16) is used to determine the tenqjera- 
ture of the isothermal section* The ten5)erature profile of 
the section II is detemdLned by using Eq, (3.12). An 
explicit finite difference sequence was used to solve 


Eq. (3.16) along wilh Bqn. (3.12). 


3.3.3 Idscussion of Besults 

The temperature profiles predicted from the above 
model are shown in Pigs. 3.10 and 3*11 > for the flow rates 
of 0.33? X 10“^ mVsec and 1.675 x 10”^ m^/sec respectively. 
It can be seen in Pig. 3.10 that the bottom temperature 



29 

remains essentially ■unaffected for about 10 hours and tliat 
the increase of the temperature of the top section is at a 
very slo-w rate. Hov;ever, for Re = 62^0, the temperature 

of top section builds up at a much faster rate. Figure 3*12 
indicates the transient rise in the surface temperature of 
the enclosed vator body for Res^iolds numbers 1250, 2500 and 
6250 . Figures 3*13 and 3«'I^ indicate the discrepancy bet- 
■ween the 'theoretically predicted and the experimentally 
measured temperature profiles. From these figures, it is 
obvious that, al'tho'ugh the one dimensional model is an appro- 
ximate one and does not allow a consideration of ■tlie flow, 
it can approximate the tenperature profile of an actual 
system quite closely, in tne ranges under consideration. 
Therefore, in several systems, the one dimensional model, 
as developed in this section, would be adequate in predict- 
ing the thermal behaviour. Tins observation is supported 
by the several studies of natural water bodies, employing 
the one dimensional model, as discussed earlier. 

3.^ RECIRCULATION IE ENCLOSURES 

The nature of mixing in bodies of water due to 
flow depends mainly on the recirculation pattern. The 
experimental study of the flovr pattern in water bodies due 
to inflow-outflow is very difficult, due to its three 
dimensicnal nature. A numerical study is carried out here. 
The two limiting cases of inviscid and creeping flow have 
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been considered initially for various geometric configura- 
tions, The results also indicate the dependence of the 
flow on the aspect ratio of the tank. In general, the 
flow is viscous and the Reynolds number is large so that 
creeping flow may not be employed. However, averaging 
in one direction, for the 2D model, particularly for low 
flow rates does bring Reynolds number low enough so that 
creeping flow may give realistic results in several cases. 
The results here are presented only as an indication of the 
nature of the flow. The actual viscous transient problem 
IS considered later, with coupling due to buoyancy. 


3,4,1 Governing Equations 

For a two dimensional, irrotational, inviscid flow 
in an enclosure, the basic governing equation, which has 
been solved by the finite difference technique, is the 
standard Laplace equations 

^ cz 0 ( 3 . 17 ) 


where is stream function. If TJ and W 


in X and Z directions respectively, then U = 
Tx 


are velocities 
ST 


and Vf = - constant at the boundaiaes and 

a constant inflow and outflow velocity gives a linear 
variation in ip at these locations. 


The basic equation for flo^'J at low Re, indica- 
tive of low velocity or high viscosity and often called 
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creeping motion, is the standard biharraonic equation 

:= 0 (3.18) 

The houndaiy conditions are obtained from the 
no-slip conditions at the walls, so that ip is constant 

3lj; 

and - 0 at the boundaries, 'diere n represents the 

direction normal to the wall. The above linear equation 
IS also the basic equation for tne hydrodynamic theory of 
lubrication. These equations are solved numerically by 
means of finite difference methods, as discussed by Gupta 
and Jaluria [ 1978 ] ^ 

3*4*2 Results and Discussions 

Figure 3*15 shows the streamlines for the inviscid 
and creeping flow, with the outflow positioned on the same 
side as the inflow, the former being at the bottom. In the 
case of inviscid flow, the streamlines are uniformly spaced, 
whereas in the case of creeping flow the stream lines are 
concentrated near the plane of the intake and outlet. 

Figure 3,16 shows the streamlines for creeping flow in a 
rectangular tank of aspect ratio 2, In the center of the 
tank, the motion of the fluid is slow. This flow pattern 
is quite different from v/hat was obtained for a square 
tank. 

Figure 3.17 sho\'/s the streainlines for a two 
flow circuit which is very common in thermal storage 
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systems. The inflow aiid o utflo\7 velocities are assumed to 
be constant in all the cases discussed above. Figure 3,l8 
sho\7s the velocity distributions at the central plane of 
the tank. The velocity W in the % plane is found to be 
maximum near 2 = 0,25 and the U component of the 

velocity IS negligible. The component w m the X plane 
is maximum at Z =0,5 and '01 is maximum at Z = 0,25 
and 0,75. 

3.? THE NONLINEAR PROBLEM 

The analysis of fluid motion within bodies of 
water due to a heated discharge is particularly complicated 
because of the buoyancy effects, which generally make it 
a mixed convection problem. In isothermal and one dimen- 
sional models, the thermal behaviour was studied without an 
analysis of the flovr, A good agreement betv/een predicted 
and experimental results was observed. This was largely 
due to the mixing caused by the flow and the buoyancy 
mechanisms. However, at lov; flow rates and constant inlet 
temperatures, such complete mixing is not expected and the 
temperature and the velocity will vary in the vertical and 
horizontal planes. A two dimensional model, assuming 
uniformity in one horizontal direction, is considered in 
this section. This model will apply v/hen complete mixing 
can be assumed in the direction normal to the inflow and 
outflow directions. This would be the case for an inflow 



33 


over the entire \n cltli of the uater body and for narrow 
water bodies, so that unifoniiity exisiis across it. In 
several cases, a tv/o cUiaensional ..'odel may be employed by 
averaging the flo\' o^/er the width of the water bodj^. The 
basic problem is, obviously, throe dimensional. But the 
complicit3'‘ involved witli solving the general problem makes 
it desirable bo analyze the flo'J m terms of simpler two 
dimensional models. In some cases, horizontal models are 
suitable and allo^r a study of the flow \/ithout considering 
the mixed convection effects, see, for instance, the work 
of Jirka, et, al, [ 1975] and Schgal and Jaluria [ 19791 • 

3 . 5» 1 Governing Equations for a Two 
Dimensional Model 

Consider a v/ater body in which the motion of the 
fluid is tvxo dimensional, as shovm in Fig. 3*7 Cd), !Ihe 
basic governing equations are those of conservation of 
mass, momentum and energy. The governing equations are 
obtained for a two dimensional flow, assuming constant 
properties and making the usual Boussinesq' approximation 
[ see, e.g,, Schlichting, 1968] as: 
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3 X 3 z 
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si = “ t (3.22) 

where u is the horizontal velocity in ^ direction and 
w IS the Vertical velocity in 2 direction. Here, the 
subscript a refers to a reference value, generally the 
ambient medium or the initial value. 

By cross-differentiating Equations (3.20) and 
(3.21) with respect to z and x respectively and subract- 
ing, the pressure terras are eliminated and the vorticity 
equation is obtained as: 
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^^^here the vorticity g is given by 
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If we define the stream function ifi by u = 
and w = " » the continuity equation is satisfied 

and it is related to vorticity by 


C = 

Therefore, the governing equations reduce to 
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The equations are nonnalized by defining charac- 
teristic time T*= ~ and d \/here is the 

reference velocity, taJcen as tlie inflow velocity. The 
following dimensionl'Bss variables are used 
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The above governing equations can tlien be \:ratten 
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where the overbars have been dropped for convenience. 

The parameters A^ , and A^ in the above equa- 
bions are given by 

= Gr/Re*^j A^ = 1/Re’ ; A^ = 1/(Re’Pr) 

where Gr = g 0 AT/-^, Re’ = d/v and Pr = v/a 

Cabelli [1976] has also considered the above 

system of equations in his analysis of flow in a solar 

energy storage tank. 
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3»5*2 Boundary Conditions 

At tile velLIs the no-slip boundary condition is 
employed. Therefore, U = 0, v/ = 0, Consequently ^ 

dilj 

is constant aaid -r— = 0, 

an 

At the inflow: 

U = 1 5 W = 0 

Therefore, the stream function was distributed linearly 
across the opening. It v;as assumed thnt the fluid entering 
the solution region was irrotational and therefore inflow 
vorticity was assigned a value zero. The vorticity value 
at the outflow was determined by using Equation (3.31 )• 

The fluid enters the solution region at a cons- 
tant temperature. The tei]: 5 )erature gradient normal to the 
boundary was equal to zero since the tank was assumed to 
be insulated at the boundaries. Ho\7ever, losses can be 
taken into account as was done in tlie case of the unstrati- 
fied and one dimensional models. 

3 , 5,3 Method of Solution and Results 

The energy and vorticity transport equations 
where both solved by alternating direction implicit (ADI) 
method of Peaceman and Rachf ord as described in Carnahan 
[ 1969 ] . The stream function equation wqs solved by using 
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the Gauss - Seidel iteration riethod. The velocity QOii 5 )onent 
in the solution domain was calculated by using second order 
centered space Oifference, The vorticity at walls was 
calculated from the streamf unction field, using the method 
described by Roache I I 976 ] , 

Finite difference methods were used for studying 
the flow pat bom. Equations ( 3 . 29 ), alongwitli (3*30 , were 
initially solved for A-j = 0 , i,e, withoat considering 

the buoyancy effects, in order bo get the flow pattern. 
Figures 3.19 and 3.20 indicate the flow patterns for 
Re’ = 100* These figures indicate the effect of geometry. 

Figure 3.19 shows the flow pattern for a square tank, where- 
as Fig, 3.20 shows that for a rectangular tank, with an 
aspect ratio of tv/Of It was found that tiie flow stabilizes 
very soon, altliough the problem could not be solved for the 
actual Reynolds numbers encountered, in the experiment, due 
to stability problems. If we reduce the mesh size, it v/as 
possible to go to higher Reynolds number, but then the 
computer time also becomes much higher. The flow pattern 
obtained is physically expected and is similar to that 
obtained earlier for the case of creeping flow. 

Figures 3,21 and 3.22 indicate the str§am line 
pattern and the isotherms for S2L.- = I.O 5 Re* = 50 

Re»2 

and T = 0 . 99 . For the same parameters, the stream lines 
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patbern and isotherms are shown in Figs. 3.23 and 3.24- for 
a time t = 2,8712, From these figures it is again evident 
that for a buoyant f 3 ow also, the flov; stabilizes very soon, 
whereas the temperature penetration is very slow, Tnis 
was also seen in the experimental study, where the tempera- 
ture field was found to var;'" ovex* a considerable period of 
time , 

The major problem in such a computational study 
IS bhe instability tliat poses a limibation on the time 
step. The coupled natures of the equations makes the 
solution very slow and, oven though the con^juter solution 
IS sometimes necessary for certain practisEal problems, 
such as power plant heat rejection systems v/here the recir- 
culation pattern is needed for the design of the system, 
the stability problem at higlier Re would often maJie the 
numerical solution a formidable task. This study was under- 
taken largely to obtain the basic procedure for determining 
the flow in a two dimensicxial system and to indicate the 
basic features of the numerical results thus obtained. 

For the experimental study carried out, it has 
already been shown that one dimensional models are adequate 
in determining the thermal behaviour. Hov/ever, interest 
does exist in the flow field in circumstances where the 
recirculation is crucial to design, as in airconditioning 
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and heat rejection systems. In such cases, a numerical 
study of bhe flov/ may be undertaken as given above to 
determine the flow pattern and the velocity levels encoun- 
tered, for steady state as well as for transient cases. 



CHAPTER If 


SUHMARY AI® CONCLUSIONS 

4-.1 SUMMARY 

In several heat rejection and energy storage 
systems, the problem of thermal discharge to an enclosed 
water body arises. The thermal effects that result are 
often of importance in the consideration of the ecological 
impact of the energy rejected on the water body. The ther- 
mal stratification that frequently arises in such cases is 
also of significance in the study of thermal energy storage. 
The present work is directed at determining the transient 
thermal effects in an enclosed water body due to thermal 
discharge. 

The study is carried out experimentally, followed 
by relevant analysis. The tnermal behaviour of the water 
body IS dependent on various variables such as inflow-outflow 
locations, inflow tenperature and flow rate, dimensions and 
geometry of the water body, environmental conditions, etc. 
This study considers several of these in detail. The effect 
of the inflow- out flow locations and of the inflow tenpera- 
ture and flow rate are considered in detail, IHie temperature 
profiles are obtained as functions of time and the growth 
and decay of thermal stratification are studied in detail. 
Based on the experimental results, analytical models are 
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developed. The complex two dimensional model, with coupling 
due to buoyancy, was studied. However, it was found that 
simpler isothermal and one-dimensional models are suitable 
for the range of flow rates examined in the present v/ork. 
Good agreement between analytical and experimental results 
was obtained, indicating the validity of these simplified 
models in the circumstances studied experimentally. For 
other cases, the two-dimensional model may be employed and 
some results are indicated. The general three dimensional 
problem is a very complex one and need be undertaken only 
in cases where the sin5)ler models are inapplicable. This 
study considers these various models and the analytical 
results thus obtained in terms of the experimental results. 

4,2 CONCLUSIONS 


The experimental results, along with various 
suitable mathematical models, have been discussed in detail. 
The results from analysis have also been compared with the 
experimental results, indicating good agreement. The main 
conclusions of bhe present investigations are disciissed 
below. 

The experimental work indicates that for the 
large flow rates considered here, Reynolds numbers ranging 
between 2000 and 62^0, complete mixing in the flow region, 
between the inflow and outflow, arises due to the flow and 
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the "buoyancy effects. For such flows, the geometric loca- 
tion of the inflow and the outflow has an important effect 
on the thermal performance of the system. If the inflow 
and the outflow are both located at the top, but on opposite 
sides, the water body is stably stratified and the level of 
stratification increases with time. If the inflow is located 
at the bottom, with the outflow at the top, the vra.ter body 
IS found to be isothermal. If the inflow is located at any 
other location with the outflov; at the top, the enclosed 
water body can be divided into two distinct regions, the 
upper region being isothermal and the lower stratified. In 
addition to these, the flow configurations commonly encoun- 
tered in various heat rejection systems, such as those in 
thermal power plants, have been studied. The effect of the 
flow variables on the outflow tomporature and on the transient 
temperature profiles have been obtained and discussed in 
terms of the physical processes involved. 

Suitable theoretical models have also been consi- 
dered for determining the thermal performance of an actual 
system. Isothermal models were found to characterize the 
transient behaviour of water bodies in which complete mixing 
arises due to flow. These were considered for constant 
inflow temperature and for constant energy input. The 
effect of energy loss to the environment on the theimal 
performance was studied. For stratified water bodies, 
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obtained wi15i specified locations of inflow and outflow, 
relevant one- dimensional models v/ere developed and the 
goveimng equations were solved numerically, Theoretically 
predicted results have been compared with the experimentally 
observed values and a good agreement was observed. 

The two dimensional vertical flow pattern has 
been obtained for the two limiting cases of inviscid and 
creeping flow. The effect of geometry on the flow has been 
investigated. The actual non-linear viscous problem has 
also been considered. It was observed that the flow pattern 
is established relatively fast and thermal diffusion conti- 
nues more gradually. The problem could be solved only for 
low Reynolds number, because of numerical instability prob- 
lems, A reduction in the mesh size makes it possible to 
go to somevdxat higher Reynolds numbers. But this involves 
large computer times and the con^jlexity increases as Re 
IS increased. Since the basic mechanisms could be studied 
at low values also, it was decided not to pursue the work 
to large values. Moreover at large values, the one dimen- 
sional models were found to be quite satisfactory, 

4,3 FUTURE WORK 

This study has considered the themal effects 
in an enclosed water body due to heated discharges at 
fairly high values of the inflow Reynolds numbers , The 
actual problem is, obviously, three dimensional. However, 



for small water bodies, conqplcte mixing may often arise 
leading to the isothermal approximation. Similarly, in 
narrow water bodies, complete mixing may exist in one or 
both horizontal dimensions, leading to two or more dimen- 
sional models. Future work must examine this question in 
detail. This would involve consideration of a vory v/ide 
range of flow rates and/or \rater bodies of various dimen- 
sions. The one-diraonsional model is a very simple one and 
has, therefore, been used for natural water bodies exten- 
sively. It IS necessary to determine its validity. The 
present work indicates some of the circumstances under 
which it may be employed. 

Another question "vdiich is of considerable signi- 
ficance pertains to the relative importance of buoyancy 
effects. The parameter Gr/Rc must be varied to determine 
how buoyancy affects the thermal behaviours. In the present 
work, this parameter is generally quite small. At largo 
values, the buoyancy- induced flov/ will dominate and the 
stratification profiles may be quite different. 

Considerable v/ork needs to be done on the under- 
standing of the flow that arises in a general three- 
dimensional problem. The two dimensional problem must also 
be studied at larger Re, With thermal stratification and 
dominating turbulent transport, the eddy viscosity and 
diffusivity will be strong functions of the depth and this 
must be incorporated in the analysis, as done for stratified 
natural water bodies. 
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